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q^' ' General Physics Institute, Russian Academy of Sciences, Moscow 1 19991, Russia 

o 

■ Key words Graphene, magnetoplasmons, collective excitations. 

PACS 73.20.Mf, 71.35.Ji 



We study optically-induced collective excitations of graphene in the presence of a strong perpendicular mag- 
netic field and a single impurity. We determine the energies and absorption strengths of these excitations, 
which become localised on the impurity. Two different types of impurity are considered i. the long-range 
^ ' Coulomb impurity, ii. a 5-function impurity located at either an A or B graphene sublattice site. Both im- 

, purity types result in some bound states appearing both above and below the magnetoplasmon continuum, 

although the effect of the short-range impurity is less pronounced. The dependence of the energies and 
oscillator strengths of the bound states on the filling factor is investigated. 
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O . 1 Introduction 
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Graphene, an atomically thick layer of graphite, has received a tremendous amount of attention since 
it was initially isolated five years ago [1]. From a theoretical perspective, this is largely due to its linear 
dispersion relation at low energies, which occurs around two inequivalent corners of the Brillouin zone, the 
K and K' valleys. Close to these points, the single particle states have a spinor character and well defined 
chirality due to a pseudospin resulting from the inequivalent A and B sublattices of the honeycomb lattice. 
There has been much discussion in the literature regarding the role of electron-electron (e-e) interactions 
in graphene, particulary in relation to disorder [1]. In undoped graphene the most striking effect of the 
long-range Coulomb e-e interaction on the material's low energy properties, is simply a renormalisation 
ON , of the Fermi velocity. This and the irrelevance of short-range e-e interactions are well established results, 

which have promoted the idea that e-e interactions are largely unimportant in graphene. However, recent 
work has portrayed it very differently, as a strongly interacting quantum liquid [2]. In the present study, 
we examine the optically-induced excitations of monolayer graphene in a strong perpendicular magnetic 
field, which become localised on a single impurity with an axially-symmetric potential. Single particle 
excitations are mixed via the Coulomb interaction resulting in a collective excitation of the whole system. 
There is little discussion in the literature of the effect of a strong magnetic field upon the e-e interactions, 
although optical properties have been suggested to be sensitive to them [3]. They are treated here beyond 
the mean field level and are central to the problem, significantly altering even the qualitative nature of the 
bound states. We explore how these states are influenced by the filling factor and the nature of the impurity, 
considering a <J-function scatterer and a Coulomb impurity. The effects of the light polarisation and the 
impurity strength are presented in Ref. [4,5]. 
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2 Theoretical Model 

Let us first recall [4] the single particle problem in the absence of an impurity potential. A single electron 
wavefunction in, e.g. the K valley (pseudospin a =ff), is *„ s ^ m (r) = (r|c* s ^ m |0) = $„f m (r)x s = 
a n (Sn4>\n\-im(' r )> <A|n| m( r )) 0, 0)xs, where the symmetric gauge A = x r is used to describe the 
perpendicular magnetic field. Here, n is an integer LL number, </> ra m(r) is a 2DEG wavefunction with 
oscillator quantum number m = 0, 1, . . ., a„ = 2^ < - SnA '^ 1 \ S n = sign(n) (with So = 0) and \s is the 
spin part corresponding to two spin projections s = j, {. The corresponding wavefunction in the K' valley 
(<r =^J.) is obtained by reversing the order of the spinor components. We consider Landau levels (LLs), 
which are split into four non-degenerate sublevels, due to Zeeman splitting (huig) an d an additional valley 
pseudospin splitting (hu> v ) with hio v <C hoj s (cf. Fig. |T}. Valley splitting occurs in high magnetic fields 
and has been seen clearly for the zeroth LL, with evidence for the n = ±1 LLs observed, but somewhat 
weaker [6]. Defining the composite indices N = {riser} and N = {nam}, the single particle energy is 



e^f = S n ^c\/\n\ + fi^sSz + ftid v cr z , where h\o c = vpy^ehB / ' c is the cyclotron energy in graphene. 

The impurity is described by an axially symmetric potential, which is given by Vc(r) = ±e 2 /ej mp |r| for 
the Coulomb impurity, where £j mp is the effective dielectric constant for the electron-impurity screening. 
Impurities with charge ±e are in the subcritical regime, so that screening effects due to the electron system 
and the media surrounding the graphene layer can be modelled by an effective charge or dielectric constant. 
Its long-range nature means that is does not scatter between the valleys and its position i.e. whether it is 
located at a lattice site or at the centre of a hexagon is largely irrelevant. In contrast the short-range impurity 
does scatter between the valleys and we consider such impurities situated at a lattice site. For an impurity 
located on an A or B sublattice site at the origin, the correction to the Hamiltonian is 



Vx(r) = V 



/ 8(r) tf(r) \ 



5(r) <5(r) 

\ / 



V B (r) =V Q 



( 

8(r) -5(r)u* 



\0 -5(t)u> <5(r) J 



(1) 



respectively, where w = e 27 ™/ 3 [7]. The Coulomb potential is diagonal in both the sublattice and valley 
indices. In contrast, the structure of Eq. (T} imposes complex selection rules, determining which transitions 
are connected via the short-range impurity. 

We consider excitations where an electron is promoted from one of the uppermost filled states A/2 
to an empty state in a higher lying LL Afi leaving behind a hole. We introduce operators of collec- 
tive excitations as Q^ 2Mz = Emi,m a =o Aa/iA^m,, (mi, m?) ^ imi d\f 2m2 , where the hole represen- 
tation, cjsi m — ► d\f and c\r — > dj^ m is used for all filled levels and the expansion coefficients satisfy 
A/ViAiMz (mi, m.2) ~ ^Af.lrnl— mi-lnal+ma- The general secondary quantised form of the Hamiltonian is 

H = ^2 \ S W S mm'eM + VtATm ) (cjv'm'^™ ~ d \f'm' d ^m) (2) 

Af \M' m,m' 

i V"^ v\ lA^m'jA^m!) t jt J /o-j 

M[,N 2 m 'l' m 2 

The first term gives the single particle energies (e/v), which are lowered by an exchange self energy cor- 
rection together with the impurity interaction. The impurity matrix elements are given by Vij^nT = 
6 SS ' J J dr 2 $j v ,(r)l^$jv(r), where i £ {A, B, C} denotes the impurity type. The second term gives the 

2 

electron-hole (e-h) interactions governed by the Coulomb potential, U (|ri — r^l) = £ | ri e _ r ,| , where e is 
an effective dielectric constant for e-e interactions in graphene. It is made up of the e-h attraction and the 
exchange repulsion, W,/™ 1 ,/ 2 !™ 2 = — W./™ 1 .^™, 2 + i^ 21 ™ 2 ^" 1 , 1 where the definition of matrix element 
U is specified elsewhere [5]. 

Copyright line will be provided by the publisher 



adp header will be provided by the publisher 



3 




Fig. 1 Diagram indicating mixed excitations for a sys- 
tem with v = 1 illuminated by a + polarised light with a 
<5-function impurity on the A sublattice. The dashed box 
encloses the excitations which are mixed for the case of a 
Coulomb impurity. 




Fig. 2 Bound states plotted against the filling fac- 
tor v for a system illuminated by a + polarised light 
with a charged donor Coulomb impurity, D + . The 
dashed lines act as a guide to the eye. 



In order to solve the problem numerically, we need to decide which transitions contribute significantly 
to the optical response of the system and to reduce the number of terms in (O to a finite number. In our 
calculations we assume all LLs with n < are filled, all LLs with n > are empty and that the sublevels 
of the zeroth LL become successively completely filled (v — 1, 2, 3, 4). Our focus is on determining the 
optically bright localised collective excitations. Only single particle transitions with no spin or pseudospin 
flips, no change of oscillator quantum number and \n\ \ — \ni | = ±1 are optically active in the two circular 
polarizations <j ± . This means that only transitions with M z = ±1 are active in ct ± . It can be shown that 
the only mixing of transitions with no spin or pseudospin flip by the e-h Coulomb interaction is via the 
exchange interaction to other transitions with no spin or pseudospin flip. We only include transitions which 
are in resonance with each other, neglecting the weak correction of higher order terms, so that for the filling 
described, only the transitions from LL n = to LL n = 1 (0 — ► 1 transitions) and —1^0 transitions 
are considered. For the Coulomb impurity there are four transitions that must be included for each filling 
factor; the relevant transitions for v = 1 are shown within the dashed box in Fig.Q] The situation for the 
^-function impurity is somewhat more complicated, since it mixes different transitions, which may have 
different M z quantum numbers. The number of relevant transitions changes with the filling factor, light 
polarisation and the sublattice, which the impurity is located on. An example of the important transitions 
for the cr + light polarisation, v = 1 and an impurity on the A sublattice is shown in Fig.Q] We diagonalise 
the Hamiltonian numerically and obtain the energies and eigenvectors of the bound states. This requires 
cutting off m at a finite value, which is justified as we seek excitations which are localised on the impurity. 



3 Results 

In Figs. |2] and [3j we show the energies of the bound states (in units E$ = (7r/2) 1 / 2 e 2 /e£b, where Ib is the 
magnetic length) plotted as a function of filling factor v = 1,2,3,4 for the case of the Coulomb impurity 
and an impurity on the A/B sublattices respectively. The Coulomb impurity is a donor D + with charge 
+e; the <5-function impurity has strength Vq = Wa 2 , where a is the lattice constant and W = 300eV. In 
both cases the system is illuminated by light in the er + polarisation. In pristine graphene all the collective 
excitations are of an extended nature and possess a quasimomentum K. They are termed magnetoplasmons 
and form a continuum of width ~ Eq contained here within the shaded area. The impurity results in some 
states becoming localised and splitting off from the magnetoplasmon continuum. For the case of the 
Coulomb impurity, we have also found quasibound states within the band, which have a high probability to 
exist on the impurity, but also long-range oscillating tails. Signatures of such states were also seen for the 
^-function impurity, but were not as pronounced. The symbol size in Figs. [2] and [3] is proportional to the 
oscillator strength (|d 2 |), with a magnification factor indicated in the legend; diamonds represent the dark 
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Fig. 3 Bound states plotted against the filling factor v for a system illuminated by a + polarised light with a 5-function 
impurity on a the A sublattice and b the B sublattice. In real systems there will be a mixture of both. 



states. Notice the differences between the graphs for the short-range impurity on the A sublattice and that 
on the B sublattice. The results are the same for even filling factors, but not for odd filling factors. This is 
because the valley splitting means that for odd filling factors the valleys are not equally filled making them 
inequivalent, which in turn introduces an inequivalence between the A and B sublattices. A significant 
difference between the impurity types is that there are fewer states localised on the Coulomb impurity 
and they evolve more smoothly as a function of the filling factor than those localised on a short-range 
scatterer. They form branches, which are indicated by the dashed lines in Fig.|2]with oscillator strengths 
which behave monotonically with increasing filling factor. The reason for this difference in behaviour is 
rooted in the greater complexity of the transitions involved for the (5-function impurity and the fact that the 
number of relevant transitions and their M z values change as a function of filling factor v. In contrast, for 
the Coulomb impurity exactly four transitions are relevant for every filling factor each with M z = 1 and it 
is only the distribution between — > 1 and — 1^0 transitions that changes. 

4 Conclusions 

We have studied the formation of optically-induced excitations bound on both a long-range (charged) and 
short-range (neutral) impurity in graphene in the presence of a strong magnetic field. The Coulomb inter- 
action between different excitations mixes them, resulting in collective excitations of the whole system. 
In a real system, the short-range scatterer can be thought of as an impurity in the graphene layer at one 
of the lattice sites, whereas the Coulomb impurity is more likely to represent an impurity in the substrate 
very close to the graphene layer. For such states to be detected by magneto-optical spectroscopy, large 
clean samples are required so that the impurity density is low enough to be modelled by our single impu- 
rity model, whilst simultaneously the number of impurities is high enough to produce a detectable optical 
signature. However, with the current high focus of research in this area, we are confident that such samples 
will soon be available. It is important to note that the bound states due to a Coulomb impurity are more 
likely to be experimentally detectable, since the <5-function impurity is less effective at producing bound 
states (they are only formed at high impurity strengths), as detailed in a separate work. 
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